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Abstract 

In this paper we study in detail different types of topological solitons which 
are possible in bilayer quantum Hall systems at filling fraction u = 1 when spin 
degrees of freedom are included. Starting from a microscopic Hamiltonian we 
derive an effective energy functional for studying such excitations. The gauge 
invariance and CP^ character of this energy fuctional and their consequences 
are examined. Then we identify permissible classes of finite energy solutions 
which are topologically non-trivial. We also numerically evaulate a represen- 
tative solution in which a pseudospin (layer degrees of freedom) bimeron in a 
given spin component is intertwined with spin-skyrmions in each layer , and 
and discuss whether it is energetically favoured as the lowest lying excitation 
in such system with some numerical results. 

I. INTRODUCTION 

Systems which permit topological excitations, i.e. where field configurations can be 
classifield by homotopy sectors characterised typically by some winding number, have been 
studied in a general sense in mathematical physics for a long time. That such interesting 
possibilities can actually arise and play a significant role in the Quantum Hall physics was 
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demonstrated in the work of of Sondhi et al [0. They showed that for example , in a 
single layer Hall liquid at filling factor = 1, the lowest energy excitations in spin, for low 
Zeeman coupling are the so-called Skyrmions and not single spin flips. These Skyrmions 
are topological excitations in the spin texture, in which the spin starts being, say, "up" 
at the origin and as you go outwards, starts tilting down in a flared manner to become 
asymptotically "down" spin at large distances. Subsequently experimental support for the 
existence of such excitations was also discovered in NMR measurements . 

Meanwhile Quantum Hall phenomena have also been studied in double layer systems , 
Q. The double well Hall plateaux at unit filling can be understood by associating with each 
electron a "pseudospin" in addition to its lowest Landau level (LLL) orbital wavefunction 
0, 0. The up and down components of this pseudospinor give the probability amplitudes 
for the electron being in the upper and lower layer respectively. The ground state of the 
u = 1 double layer system, known to be a quantum Hall state with a Hall conductivity 
plateaux is a pseudospin ferromagnet, with the pseudospin aligned in the x-direction. This 
is a very remarkable phenomenon in that it amounts to interlayer coherence between the 
electrons in the two layers. This pseudospin degree of freedom is in addition to physical 
spin. To start with, in analysing double layer phenomena, the spin degrees of freedom are 
suppresed for simplicity. Even then one can still consider excitations in the pseudospin. 
Inspired by the presence of Skyrmions in spin, people have also considered the possibility 
of topological excitations in pseudospin. Such pseudospin textures called "merons" and 
"bimerons" have been suggested as possible low lying excitations of double layer systems 0. 
The homotopy group vr2[S'2] and its winding number are identical for spin and pseudospin 
since mathematically pseudospin is identical to spin, both being SU(2) spinor fields on 
a plane. The change in terminology from Skyrmions to bimerons does not indicate any 
topological difference between the two in going from spin to pseudospin excitations but only 
differences in their detailed profiles. This difference in turn happens because of the difference 
in the energetics of spin and pseudospin and correspondingly , their asymptotic direction. 
Meron excitations , if present in double layers, can give rise to a Kosterlitz-Thouless (KT) 
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1^ transition which may be enable them also to be experimentally observable. 

Clearly there are prospects of even more esoteric excitations when both spin and pseu- 
dospin degrees of freedom are considered simultaneously. That is the theme of this present 
work. We will continue to study the unit filling factor {u = 1) case. There has already 
been some discussion of the combined spin-pseudospin z/ = 1 double layer system [^], 0. 
Our work discusses diferent aspects of the problem than these studies. We analyse in sub- 
stantive detail intertwined spin-pseudospin topological excitations of this system . When 
both spin and pseudospin are active degrees of freedom , these are together described by 
a 4-component object. This 4-component object has been referred to as a CPs spinor in 
the hterature 0. That is correct, but needs to be justified. A theory does not become a 
CPs theory just because its field is a normalised 4-component object. The system must 
enjoy a U(l) gauge invariance, which is what makes the spinors span a projective space, 
implied in the acronym CP. Without that gauge invariance the beautiful results on CPn in 



the literature |^, [|TT| cannot be borrowed. 

So we begin in sec. II by showing, starting from the basic microscopic theory of the 
V = 1 system that in the effective LLL theory for the spin-pseudospin texture such gauge 
invariance is_ there. This is a straightforward derivation following the procedure developed 
by Moon et al In fact we find that in the limit where the layer-separation d vanishes, 
the Coulomb interaction energy is precisely the protoype CP^ Euclidean action used in the 
pioneering papers on that topic |jTO| , for which exact topological solutions are known in terms 



of analytic functions. Of course, when (i 7^ 0, the energy functional is more complicated and 
these analytical solutions do not hold. But the theory is still a CP3 theory, and the homotopy 
classification of the solutions still holds. Only the solutions themselves have to be calculated 
numerically. 

A topologically non-trivial CP3 solution will generally involve an intertwined texture in 
the physical spin in each layer as well as in the pseudospin of each physical spin projection. 
One can ask whether such solutions can be legitimately interpreted as containing, as subsys- 
tems, spin-Skyrmions in either or both of the layers, possibly intertwined with a pseudospin 



meron or bimeron . If so, then such possibihties of containing several topological entities 
as subsystems has to be made compatible with the fact that any finite energy CP3 texture 
carries altogether only one topological winding number. We study all these questions in 
sec. Ill and find that there are certain restrictions on the types of solutions permitted. We 
show that the individual layers of a double layer system cannot accomodate all possible spin 
structures one may find in a pair of unrelated single layers. The spin winding numbers in 
the two layers are related to one another and to the pseudospin winding number. 

Consistent with these restrictions, we then pick in sec IV a representative ansatz which 
can be viewed as a spin-Skyrmion intertwined with a pseudospin bimeron. We then numeri- 
cally evaluate such a solution by solving the coupled nonlinear partial differential equations 
that arise from extremising the texture energy functional. In earlier work |T^, |13[ we had 



studied in some detail both meron and bimeron excitations in pseudospin for double layer 
systems, with the spin degree of freedom suppresed. The present calculation is a more com- 
plicated version with CP3 spinors, but is done by similar numerical techniques. We present 
the spin and pseudospin profiles of our intertwined solutions for different values of interlayer 
separation. 

We also estimate the interaction energy of these solutions for some typical sets of values of 
system parameters. We discuss the dependence of this energy on the separation between the 
two meron centers. We find , as expected, that if only the gradient and capacitance energies 
are considered , their minimisation will drive the textures towards zero size. Therefore we also 
calculate the topological charge dependent Coulomb energy of our solutions which, being 
repulsive , should drive the merons farther apart, offsetting the above tendency towards zero 
size. Then we extremise the total energy so obtained and find that it does show a minimum 
at some optimal meron separation, for each value of layer separation. 

We also find that these energies are approximately of the same order as those of purely 
spin Skyrmions of the single layer system. We make qualitative speculations on whether 
or not our spin-pseudospin intertwined solitons can be energetically favoured over solitons 
purely in spin or pseadospin, or over simple spin-flips. 



II. TEXTURE ENERGY AND ITS GAUGE INVARIANCE 



In a double layer quantum Hall system with both spin and pseudospin degrees of freedom 
present, an electron will carry, apart from its coordinate wavefunction (pxi''^, a 4-component 
normalised spinor whose components in general may vary with the orbital quantum number 
X. For any given X, this spinor can be denoted by 

/ai(X)\ 
a2{X) 



aa{X) 



(2.1) 



asiX) 
\a,iX)J 

where the spin-pseudospin index a = 1,2,3,4 corresponds to amplitudes that the electron 
is in the upper-layer up-spin, upper-layer down- spin, lower layer up-spin and lower -layer 
down-spin states respectively. It will henceforth be understood that the spinor is normalised, 
i-e.I](T I (^a{X) p= 1 for each X In the literature, this a^j has sometimes been referred to as a 
CPs spinor (see for instance Ezawa Q]). That is correct, but requires a little justification. In 
a CP3 theory, the spinor must not only be normalised 4-component object, but be defined 
only modulo a local gauge transformation common to all four components. This in turn 
requires that the Euclidean action or static energy functional of the spinor field enjoy a 
corresponding gauge invariance. In this section we will verify all this. We will also see that 
the nature of the gauge symmetry is different for a double layer system than for a pair of 
isolated single layers. This, as we shall see, has the important consequence of prohibiting 
certain topological spin excitations in the double layer system which would have been present 
in the individual layers had they been far apart. In this way, along with establishing the 
CPs nature of the system we will also identify permissible types of excitations where the 
spin and pseudospin are nontrivially intertwined, some of which we numerically evaluate in 
later sections. 

Let us start by deriving the energy functional of any spin-pseudospin texture from the 
microscopic Hamiltonian. This is just a straightforward generalisation of the procedure 
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already in the literature for the simpler case of a spinless bilayer problem ||^. Therefore we 
need to present only the essential equations needed for completeness and understandability. 
We take the microscopic Hamiltonian to be 

H = Hk + Hi + Hc (2.2) 

Here 

Hk = ^Y^j dfijt (2.3) 

is the kinetic energy in the presence of the magnetic field. We will be working at z/ = 1 in 
the lowest Landau level (LLL) approximation. Corespondingly, the operator ip(j{r) is the 
LLL-projected electron field operator expanded in terms of Lowest Landau Level orbitals as 

N 

^.(r1 = E C'rx (2.4) 

x=i 

with 0x(^ being a LLL orbital, say, in the Landau gauge with X as its guiding center. 

The second term in the Hamiltonian is the one body term representing the Zeeman and 
interlayer tunnelling energies. 

H, = Y.I dr^ii^ i~9^^ - (2.5) 

(7,5 

where and are spin and pseudospin matrices suitably generalised as 4 x 4 matrices on 
the outer product space of spin and pseudospin. 

The third term in the Hamiltonian is the Coulomb term : 

1 ^ r 

= 2 E y ^^^lt^^l<(n)V'I,(rl)^'^^'^^(ri-rl)^.,(rl)^.,(fi) (2.6) 

o"l,a'2 = l 

In the above , the Coulomb potential depends on whether the particles are in the 

same layer or different layers 

ycricr2 _ ^ ai , a2 in same layer, 

e ri2 

e2 

^0-10-2 _ _ — ____ ^ different layers (2.7) 



where d is the interlayer distance. To obtain the energy of an arbitrary spin-pseudospin 
texture, we adopt the strategy followed in the work of Moon et al 0. We first consider the 
ansatz state 



where | 0) is the vacuum (no electron) state , X stands for Landau gauge orbitals and a^-^X) 
is an orbital dependent 4-spinor as in eq.(2.1). In the high-B limit each Landau gauge orbital 
density is uniform along the y-axis with support on a thin line localised around some value 
of X. Further these states are closely spaced along the x-direction. Using this feature, we 
will later on replace the orbital label X by the x-coordinate itself. In that case the above 
texture a„{X) depends only on the x-coordinate and not on y, and therefore carries zero 
topological number density (see eq.( 3.14) below). Nevertheless we will use this ansatz to 
calculate its energy functional, and then later that energy functional to the more general and 
topologically non-trivial textures by invoking isotropy of the system in the x-y plane. This 
was the exactly the strategy used in ref( |^). We will calculate the energy functional of the 
spin-pseudospin texture (2.8) by taking the mean value of the second quantised Hamiltonian 
in that state . 

At unit filling z/ = 1, and in the space of LLL orbitals the kinetic term Hk is just a 
constant equal to {N/2)huj ,the energy of the filled LLL band. This constant will henceforth 
be neglected. 

The Zeeman and tunnelling one-body energies yield 



The expectation value of the Coulomb interaction Hamiltonian can be convenientUy 
written in terms of the following spinorial bilinears for the upper (u) and lower (1) layers 



= n E ci^<^Ax)] 1 0) 



(2.8) 



X cr 




(2.9) 



F„(X) = |ai(X)p + |a2(X)p 



(2.10) 
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F,(X) = |a3(X)p+|a4(X)|^ 

i=l,2 



(2.11) 
(2.12) 

(2.13) 



i=3A 



On inserting He from ( p.6|) and the state \1/ from (|2.8|) straightforward algebra then gives 
us the Coulomb energy in terms of the spinors Qa- 



EcMX)] = {^\Hc\^) = {Hc)d^rect - {He) 



exchange 



(2.14) 



with 



(He) 



direct 



2 ^ 

Xl,X2 



+ (D^ - D'){Fu{Xi)FiiX2) + Fi{X,)Fu{X2)) 



(2.15) 



and 



(He) 



exchange 



E 



(2.16) 



Here 



D''''^ (X2 - Xi) - Vxf^x2,Xi,X2 

rps,dfY' V \ -[/s,d 



(2.17) 



with 



ysd 



dfidr2V'^\n - r2) X <p*x^{riWx^{r2)(t)x,{n)<PxAr2) 



(2.18) 



These direct and exchange Coulomb interaction matrix elements D*''^ and E'^''^ between two 
electrons in LLL orbitals Xi and X2, in the same(s) or different (d) layers, are exactly the 
same as were used in the spinless double layer problem by Moon et ai 0. However, the 



inclusion of the physical spin degrees of freedom is reflected in the energy expressions in 2.15 



and 2.16 , which involve all four components of the spin-pseudospin multiplet a^. 



Adding the contributions in eq. (|2.14| ) and eq.( p.9|) we get the total energy expectation 
value 
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E[a^{X)] = E,[a,{X)] + ^cK(X)] (2.19) 



In the Hartre-Fock approximation, this energy expectation value i?[ao-(X)] in eq.( p.l4[ ) will 
be minimised to get the ground state and excited state spin-pseudospin textures. 

But, let us first examine the gauge invariance of the energy functional E[a„{Xy\. Consider 
the transformation 

a„{X) e^^-Wa,(X) for a = 1,2 

a,{X) e*^'Wa,(X) for a = 3, 4 (2.20) 

Notice that we have used different phases A„(X) and A;(X) for the upper and lower layer 
components respectively. This is a f/(l) x f/(l) transformation. These phases can also vary 
with the orbital index X. [ Note : X is not the space coordinate. But, following accepted 
approximations (see ref 0) eventually the sum over the orbital index X will be converted 
into an integral over space coordinate, invoking the fact that for large magnetic fields, each 
LLL orbital wavefunction is highly localised. Hence the above X dependent transformation 
corresponds to spatially local gauge transformations.] 

Under these local ?7(1) x f/(l) transformations, the one-body Zeeman energy in ( |2.9D and 
the direct part of the Coulomb energy ( |2.15D are trivially invariant since they involve only 



the squared- modulus of ao-(X) . So is the first part (proportional to E'^ ) of the exchange 
Coulomb energy ( p.l6|) . But the tunnelling energy in (|2.9| ) and the second piece of the 



exchange energy (|2.16| ), which involves exchange Coulomb interaction E'^ between different 
layers are invariant only if 

A„(X) = A^X) = A(X) (2.21) 

Thus the full energy of the double layer system enjoys only a U(l) subgroup of U{\) x U(l) 
defined in ( |2.20| ) - a subgroup where all 4 components of a^j are transformed by the same 



phase. This is the U(l) gauge invariance modulo which our CP3 spinors are defined. 

Consider, however, what would happen if we had very widely separated v = 1 layers 
(the separation d 00). Then each can have its own 2-component spin texture described 
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by a CPi system (equivalent to a non-linear 0(3) sigma model) with its own U(l) gauge 
symmetry (see |^ and |11[). The well separated pair of layers should enjoy U{1) x U{1) 
gauge symmetry. Our derivation shows a similar effect. When d ^ oo, both the tunnelling 
parameter t in ( p.9| ) and the interlayer Coulomb potential v'^ involved in ( |2.16| ) would vanish 



and the full U{1) x f/(l) gauge invariance would indeed be restored. We will see later that 
this reduced gauge symmetry of a double-layer system at finite separation has consequences 
in terms of what types of finite energy excitations are permitted in it as compared to a pair 
of isolated single layers. 

III. GRADIENT EXPANSION AND THE CP-x FIELD THEORY 



To rewrite the energy expression p.l9| in a continuum field theory language, we proceed 
following Moon et al ^ and convert sums over the LLL label X into an integral over space. 
Clearly the one body energy ( p.9| ), which involves only a single sum over the index X, will 
become a local term, i.e. a spatial integral over the 1-body energy density. But the Coulomb 
term ( |2.14| ) containing a double sum over Xi and X2 will become a non-local term involving 



a double integral over some coordinates Xi and X2- For long wavelength excitations one then 
makes the usual gradient expansion . Expand the spinor for X2 as 

a.(X2) = a,(Xi) + ((X2-Xi)-^a^(Xi) + (3.1) 

Up till now we found the energy of textures which were y- independent. 

Now we will invoke the isotropy of the basic system in the x-y plane and generalise this 
expression for arbitrary textures by This is done by making the replacement 

and by replacing x-derivatives by gradients. Insert the above expansion ( p.l| ) into the 
Coulomb energy expressions (2.13 - 2.16). Keep terms only upto order dx_^ and replace 
the sum J2xi by an integral over space as indicated. (These steps are given in the Moon et 
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al work for the simpler spinless double layer case). The result, for our problem , is the 
following local expression for the total energy ( |2.19| ), with overall constants subtracted out : 
1 



dr 



g I |ai|^ - |a2p + los^ 



|a4|^^ — t(^aial + a2al + h.c. ^ 



+ /3^l dv{Fu{r) - Fi{r)) 
+ 2p' I dr 



i=l,4 i=l,4 



+ (p*^ - p') / dr 



(3.3) 



I 2 3*V72/ 3 2*\ I 2 4*^72/ 4 2*\ , u 

+ a a v [a a ) + a a V (a a ) + h.c. 
where the constants appearing above are defined by 
Pm = I Y {{E\X2- X,) ~ E\X2- X,)) - {D^X^- X,) - D^{X2- X,))] (3.4) 



1 X - (X2 — Xi] 

2 ^ 2 

^ (X2-X1) ^ 



p 



1 ^ {2kz^E\x,-x,) 



(3.5) 



(3.6) 



These constants are again the same as given by Moon et al in the spinless double layer 
problem. The term involving (3 represents the "capacitance energy" of the double layer 
system. It is proportional to the square of -Fu(r) — F/(r), which gives the difference in charge 
density between the two layers. The constants and p'^ represent spin-pseudospin stiffness 
coming from intralayer and interlayer Coulomb interaction respectively. 

This energy functional |3.3| will act as the effective classical Hamiltonian to be minimised 
to find different textured solutions. The ground state will correspond to a spatially uniform 
texture, and so can be obtained by minimising the gradient-free terms in p.3| . This is acheived 
by the spinor a^iX) = -^(0,1,0,1). The one-body Zeeman and tunnelling energies are 
clearly minimised by this choice since the spin is polarised "down" in both layers and the 
psuedospin is along the x-direction, i.e. a layer-symmetric state. This choice also minimises 
the capacitance energy since it has equal occupancy in the two layers Fu = Fi = ^. 
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Moving on to excited states with non-trivial textures, these are obtained by extremising 
the full energy functional (|3.3|). Note that (|3.3| ) including its gradient terms is still gauge 



invariant under the local U(l) transformation mentioned earlier, 

a.(X) e^^(^)a.(X) (3.7) 

so that this is still a CP3 theory. In fact the term proportional to the isotropic spin- 
pseudospin stiffness , namely , 

EcP = 2p' I dv[ {d^d*{r)d^a\r)) + ( ^ d*{r)d^a\r)^] (3.8) 

1=1,4 j=l,4 

is the Euclidean action for the prototype minimal CP3 theory ||10[. Indeed, in the limit 



where the layer separation d is zero, this Ecp will be the only surviving term from the 
Coulomb energy in p.3| since the interlayer and interlayer Coulomb potentials will become 
equal ( (f* = v'^)) and hence both (3 and p'^ — p'' will vanish. 

The properties of this prototype CP3 system and its topological solitons are well known 
pi]| , |TT| . Let us briefly recall those salient features which will be of relevance to us. Define 



a gauge field as follows. 

A^ = iY,[ald^a^] (3.9) 
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Clearly under the gauge transformation , 

^A^ - d^A (3.10) 

The energy Ecp can then be written in a manifestly gauge invariant manner as 

4 2 

Ecp ^ 2p^ [ dr[J2 \ I' (S-H) 

(j=i ^=1 

where = df^ + iA^ is the covariant derivative of the U(l) gauge transformation. Then 
any finite energy field must obey, as r ^ 00, the boundary condition 

D^a^ = {d^ + iA^)a^ = (3.12) 

Since A^ is independent of the spinor index cr, this implies ( see [0) that as r — > 00, 
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(3.13) 



where is some constant spinor . The important point is that all four components of 
have the same asymptotic phase (f) which may depend on the spatial angle 9. The underlying 
reason is that the system is invariant under the same single U(l) gauge transformation |3] 



acting on all the four components of a^. Finally, the phase function e 



as r 



OO IS 



a mapping of one circle (spatial infinity) into another (the U(l) group manifold), and can 
therefore be divided into homotopy classes characterised by a winding number 



Q = -77- / d^r tf,y{D^a„y {D^,a^) 
/vr J 



(3.14) 

For more details supporting these results see ref [|TT| . Exact soliton solutions for the minimal 
CP3 system also known analytically in terms of analytic functions. Those will not however 
hold for our full system p.3| which has to be used when the layer separation d ^ Q. The solu- 
tions will have to be obtained numerically by using appropriate ansatz. But the boundary 



condition |3.13| and the winding number classification will still hold. They can be used to 
decide what forms of intertwined spin-pseudospin solitons are permitted in double layers. 

An important consequence of the common phase boundary condition ^.13| is that certain 
spin textures one can imagine having for two separate single layers are not permissible in 
the double layer system. Consider a single layer aX v = 1 carrying a Skyrmion with winding 
number n. This is a finite energy configuration which can be described by a 2-component 
spinor, say, 

^ A(r) ^ 



V 



/(r)e*"^ 

obeying boundary conditions as r — > 00 given by 



A(r) 




1 



and as r — > 
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A(r) 
fir) 



1 




(3.15) 



One can have two such layers, widely separated, with two different spin-winding numbers 
n and m respectively. Nothing prohibits this . However suppose the two layers are part of 
a u = 1 double layer system at finite d , and are described by a CP3 4-spinor 

/ Ai(r) \ 
/i(r)e^'"^ 
A2(r) 
V/2(r)e^"V 



1 
71 



(3.16) 



This would violate the condition |3.13| since asymptotically the second and fourth components 
would have different phase functions. Such a texture is forbidden as per our analysis and 



indeed if one calculates its energy by inserting it in |3.3| one will find the energy diverging 
logarithmically. The divergence comes from the angular derivative of -^Oq contained in the 
Laplacians in |3.3| . That yields a contribution to the energy density proportional to 

(n + m)^ 



O 9 

n + m 



(3.17) 



as r 00 which will lead to a logarithmic divergence unless n = m. At the theoretical level 
the reason for this can be traced to the reduction of gauge symmetry discussed earlier, from 
f/(l) X f/(l) to U{1) when two layers are together. 

Keeping in mind this constraint of equal spin-winding numbers in each layer, let us il- 
lustrate non-trivially intertwined spin-pseudospin configurations with the following example 
that is allowed : 

/ Ai \ 

z-b 

A2 
\z + bj 

Here Ai,2 and b are non-zero constants while z is the complex coordinate on the plane. 
A = (A^+ \ z — b l"^ +A2+ I 2; + 6 P ig ^]-^g normalisation factor. Asymptotically, the 



A 



(3.18) 
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first and third components of |3.18| botli beliave as -^e*^ wliile tlie otlier two components 
vanisli. Tliis is tfierefore a permitted (energetically finite) CP3 configuration with winding 
number Q = 1. 

One can see that this example is so designed that within each layer the spin texture 
looks like that of a single Skyrmion, while at the same time it is also a "bimeron" in the 
"psuedospin of the down-spin component" (contained in the second and fourth components 
of the 4-spinor of p.l8| . See ref ( |jl3|) for more on bimerons.) But, we should remember that 
the the upper and lower layers are not separately normalised in the example |3.18| . As r varies 
so does the relative charge density in the two layers. Thus the spin vector in the upper (or 
lower) layer in |3.18| will not be a unit vector at every point unless it is is locally renormalised 
by the charge density of that layer at that point. Similarly, while the pseudospin of the 
down-spin component in the example |3.18| forms a bimeron, this pseudospin will also be a 
unit vector at each r only after being renormalised by the down-spin density, which varies 
from point to point. Such renormalisation can be achieved by writing any general CP3 
4-spinor (|2.1|) in terms of spin and pseudospin polar angles. 



/ cos f cos 1^ \ 



cos f sin 'fe"^ 



sin I cos |-e*^ 



Vsinf sin|e*(^+'^') 



(3.19) 



where the angles O^.i, and (j)u,i are the polar angles of the spin in the upper (lower) layer 
while a and (3 are the polar angles of the pseudospin, each of which is a the function of 
the coordinate r*. (Recall that CP3 spinor has 6 real gauge invariant degrees of freedom. 
) Suppose we tentatively define, using these polar angles the familiar expression for the 
spin-Skyrmion number in each layer by 



nu,i = ^ I ^ re^^d^{cos6u,i) d^y(pu,i) 



(3.20) 



One can then verify that the configuration |3.18| indeed yields unit spin-winding numbers 
= 1 in each layer. 
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Similarly , to get the pseudospin winding number one uses an alternate parametrisation 
of the same 4-spinor. 

/ cos ^ COS Y \ 



cos ^ sin ^e"^" 



sin ^ cos ^e*^T 



(3.21) 



Vsin^sin^e^(^i+'^=) / 

Then the pseudospin winding number for the down spin component, for example, can be 
written as 



npsii) 



1 

A-n 



(3.22) 



Again, the example |3.18| happens to yield ^ps(i) = 1 in addition to, as we have seen, riu^i = 1 . 
Thus the example |3.18| illustrates an intertwined spin- pseudospin topological configuration, 
containing the spin texture of a Skyrmion in each layer and the pseudospin texture of a 
bimeron in the down-spin component. 

One should however be cautioned that there is only one true topological charge Q in 
the full CP3 theory, given in eq. ( ^.14[ ). Although the above example ( ^.18| ) contains the 
texture of two Skyrmions and a bimeron, its CP3 topological index Q obtained by inserting 
it into eq. ( glil) will come out to be not 3 , but unity. The separate sub-charges for spin 
and pseudospin defined in ( |3.20|) and (|3.22|) in general do not have the same sanctity as they 
would have had for Skyrmions in a single layer or bimeron in a spinless problem. Although 
in the above example these separate spin and pseudospin winding numbers turn out to be 
integers, in general they need not be integers , or more importantly , be conserved in time. 
They are not protected by homotopy considerations in our full 4-component theory. The 
angles 6u^i etc used in |3.19| cannot always be obtained from the original components Oo- of the 



4-spinor p.l| , since they are not defined at those singular points where a = vr, respectively. 
A similar remark holds for the other angles used above. The numbers nu,i and 
change in time due to leakages through such singular points . It is however interesting to 
note that the exact CP3 winding number can be rewritten in expanded form using the angles 
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defined in |3.19| into parts that can be attributed to winding of spin and pseudospins. Tliis 



also brings out tlie intertwining of spin - psedospin texture . We liave, 

Q = ^J dre^''(^d^{cosa)[{^{l - cos^„)9,0„ - ^(1 - cos^z)5,0; - 8,(3] 

-F^{f)d^icos9u)d,(j)u - Fiir)d^{cosei)d,(Pi]^ (3.23) 

where -Fu(r) = (1/2) (1 + cos a{r)) and Fi{f) = (1/2) (1 — cosa(r)) are respectively the same 
quantities as in (|2.1CI| , |2rTT| ) and denote the number density in the top and bottom layers. 
For the spinless (spins fully frozen) case one can set 6u,0i,(f)uO'nd(j)i to be constants. Then 
one will recover the pseudospin topological charge formula 



Hps = / dre'"'d^{cosa)d^p. (3.24) 



Similarly for a single layer (say the upper layer) case one can set a = and recover the spin 
winding number formula 

n„ = 1 dre''''d^{coseu)d,^u. (3.25) 
In the general, where both spin and pseudospin havesome intertwining texture, the full 



topological charge will receive contributions from the windings of all these, as given in |3.23 



Finally, the very simple example p. 18 not only illustrates a nontrivial intertwining texture 



it is also an exact solution of the prototype CP^ theory p]8| since its components are 



analytic functions (see ref |11|). But, for our full theory in the presence of a non-zero layer 



separation and with Zeeman and tunnelling energies , classical solutions minimise the full 



energy functional p.3| , have to be obtained numerically by solving the non-linear coupled 
partial differential equations that the minimisation conditions yield. The simple analytical 
example |3.18 however will guide us in setting up the desired ansatz for the numerical solution 



with appropriate boundary conditions so that intertwined textures which are nontrivially 
wound in both spin and pseudospin can be obtained for our full theory. An illustrative 
family of such solutions is obtained in the next section. 
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IV. FIELD EQUATIONS AND THEIR SOLUTIONS 



Classical solutions that minimise the full energy functional 3^ have to be obtained 
numerically. To do this we use the parametrisation of the spinor components of the form 



/ cos f COS 1^ \ 



COS -^sm-^e 



\ sin ^e**^ 



i4>u 



(4.1) 



A 



(4.2) 



One can see that this is a sub-family of the general case in eq ( p.l9| ) where for simplicity 
we have set di equal to vr and f3 = 0. 

We will look for numerical solutions which would have corresponded ,if the energy had 
been of the simple prototype functional ( ^.181 ), to its exact analytic solution 

/ A \ 
z — b 


\z + bj 

This configuration represents a spin skyrmion in the upper layer intertwined with a 
bimeron in the " pseudopsin of the downspin component" . It does not have any non-trivial 
winding in the real spin of the lower layer (though the fourth component in the spinor varies 
over the coordinate space the spin will be always down). To ensure that our numerical 
solutions have the same topological properties as well as similar profiles as the prototype 
spinor ( ^.2| ) we impose the same boundary conditions on the components of the former as 
obtained in the latter, both asymptotically and at the meron centers centres x = ±b. 

In terms of the ansatz |4.1| the local energy functional ( |3.3| ) takes on the form 

a 



Ec = Pm dfcos^ a + 2p'^ / df 



i((Va)^ + cos^^(W^ 
+ -(^(1 + cosa)(l - cos6'„)(V0„)^ + 2(1 - cosa)(V0i 
- ^(^(1 + cosa)2(l - cosOuYiVcPuf + 4(1 - cos a)' {V (Pi)' 
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+ 4(1 - cos^ a) (1 - cos eu)V<pu ■ ^<p}j 
+ {p'-p') j dr\]^{{V cos af-{Vaf) 

+ 2{l-cosdu)V<pu-V<P?j) (4.3) 
This energy functional has to be minimised with respect to all the angle fields in the 



anstatz. As we did in our earlier work on the spinless problem [13|, here too we will use the 



bipolar co-ordinate system |T6[ to describe the spatial plane. 



ln\z — a\ — ln\z + a\ ; = arg{z — a) — arg{z + a) (4.4) 



We have already elaborated in the advantages of this co-ordinate system when one 



has to impose the bimeron -type boundary conditions. However here the advantages of 
introducing such an unfamiliar co-ordinate system is not as much as in the simple spinless 
bilayer problem of [|13| because the ansatz here is not symmetric between the two layers. 
Consequently unlike the spinless case cos a is no more antisymmetric about rj = axis. All 
these features along with the fact that the energy minimisation unavoidably requires solving 
coupled non-linear partial diffential equations ( p.d.e.) render the numerical exercise much 
more complicated here . What we have done under these circumstances is the following. 

We have solved the field equations numerically for the case where just the capacitance 
term is added to the minimal CPs energy . /^From our earlier calculations we know that 
this term is going to change the solutions considerably. The terms in each equation with the 
coefficient (p* — p'), which accounts for the anisotropy in the exchange energy is not included 
in the process of numerical integration. As a justification of such simplification we can say 
that the anisotropic terms which involve higher order gradients of the spin pseudospin field 
will have less pronounced effect compared to the capacitance term on the solutions. This 



has been graphically shown in Fig. 1 and 2 of our earlier work [jl3[. Even after this drastic 
simplification we are still left with solving four coupled non-linear p.d.e's . For example, the 
equation which is obtained by extremising the energy w.r.t cos a is 
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5E, 



c 



5 cos a J ps=pi 



'^•PmQs COS a + p'* 



- 2 



cosa;(V cos a) 



- 2 



cos a 
1 — cos^ a 



^1 — cos^ a)2 
l-cos^^J(V0„)^-2(V0O' 

2/V7 i \2 



-(1 + cosa)(l - cos6'„) (V0U 
1^1 — cosa)(V0i)^ + cosa(l — cos 9u)(V4>u ■ ^(pi 



where 



(4.5) 



(4.6) 



(coshr] — cos 0)2 

is the Jacobian of this coordinate transformation and all gradient operators are defined in 
the bipolar-coordinate system |T^ . Similarly we wil have three more equations obtained by 
extremising the energy functional with respect to 6u, (pu, 4>i and then writing the resulting 
equations in bipolar co-ordinates. We will not display them here. 



A. Numerical Procedure 



The numerical procedure is almost same as that in Ref. [13|. Here also one can see that 



the Jacobian factor Qg in the first term of the eq. |4.5| is singular at the point (r^ = 0, = 0). 
The behaviour of cos a near this point is also going to be same as the behaviour of in 
ref. ( [|13|). The major difference compared to the earlier problem however comes from the 
fact that it is no longer sufficient to find out the solutions in one quadrant and get the rest 
from symmetry considerations. We have to solve this problem on both side of the t] = axis 
since our starting ansatz solution is not completely antisymmetric around rj = 0. During 
the numerical work one also has to be careful about the different branches of the angles 
(pu,!- As one needs to integrate the equations on the both sides of the r] = axis the size 
of the mesh on which we have to discretise the field equations becomes larger compared to 



the earlier case of spin-frozen double layer problem |T^. Also here we have to solve four 
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coupled p.d.e's simultaneously. This simultaneous increase in the number of lattice points 
as well as independent fields demands that we have to invert a huge determinant in the 



Newton-Raphson procedure [|17| while improving over the initial guess solution. This forces 
us to increase the lattice constants of the mesh slightly compared to what we have done in 
our earlier work ||TB[. However we have checked that the error introduced in this way is not 



very high. In the next subsection we shall present our results along with the discussion. 

B. Results and Discussion 

Our solutions of eq.( [4.5|) along with the other three field equations yield the spatial 
behaviour of the CP3 fields parametrised in term of the angles a , /3 fiu.i and (\)u,i- The 



calculations are done iteratively. We start with the simple analytical spinor [4.2| which is 
the exact solution when the capacitance and anisotropy terms in the energy ( [4.6| )are absent. 
Then the capacitance term in the equation is introduced in small steps and the corresponding 
solution obtained numerically. We have performed several calculations each starting from 
different initial values of the constants A and b . The constant A represents the starting value 
of the first component of the spinor in the iteration process. It stands for the spin-Skyrmion 
size in the CP3 limit, but when subsequent iterations are performed in the presence of other 
energy terms, it is replaced by a space dependent solution. But the parameter h is fixed for 
a given calculational run. It represents the meron separation and enters into the equation 
( [4.5|) explicitly through the first (capacitance) term. While we do calculations for different 
values of 6, the optimal value of h will have to be obtained by minimising the energy with 
repsect to it. We will return to this point later. 

We present below the salient features of our numerical results. The major feature we want 
our numerical solution to have is the intertwining of the spin skyrmion with the pseudopsin 
bimeron . We would also like to show the leakage of electrons of either spin from one layer 
to another as we move in space, as a fallout of this intertwining. To show this , we have 
plotted both cos a as well as cosaj as a function of r] for a set of values of the angle in 
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Fig.l and Fig. 2 respectively. These solutions correspond to layer separation d = 0.6/ , 
bimeron center-separation b = 2.5/ and starts from an initial value of A = 1 in the starting 
trial solution |4.2| . The sequence of curves shown correspond to (p equal to O.OQvr, O.SGvr, 
O.GStt, and O.OOvr respectively with the outermost one belonging to (p equal to O.OOtt. As we 
have discussed earlier [|T3|, Note from the definition of the bipolar coordintaes that spatial 




infinity in x-y plane corresponds to f] and both equal to zero. As we approach this point 
in the 0, plane, the solution should damp exponentially as exp( , ) where 

Mi (4.7) 

PA 

. Correspondingly we see in Fig. 1 and 2 that the low curves rise very slowly as r] increases 
away from zero. 

The interesting point to note about these solutions is that in Fig. 1 cos a approaches dif- 
ferent (absolute) asymptotic values as r] approaches ±oo. (These are respectively the centers 
of the two merons that form the bimeron. Although computational limitations allow us to 
go only upto values Of r/ = ±3, it is clear from the figure that asymptotic behavior has been 



obtained ). This asymptotic behaviour is extracted directly from the analytic ansatz ^]2 
and implies the leakage from the pseudospin to spin. It is useful to remember at this point 
that we have a bimeron only in the "pseudospin of the down-spin componenet" whereas 
cos a represent the z-component of the total pseudospin. This is realised in Fig. 2. Here 
cos a i{down — spin) represent the z-componenet of the "pseudospin of the down-spin com- 
ponent". It is completely antisymmetric about t] = and approaches ±1 as 1] approaches 
±oo. This behaviour is same as the behaviour of in the spinless bilayer case. This is how 
we can extract from our results the pure bimeron by suitably partitioning the pseudospin 
into different spin components. 

Since bipolar co-ordinates are not very familiar we have given an alternate representation 
of the above results through a vector-plot in the physical x-y space in Fig. 3 and 4. The 
values of the parameters in these figures are same as those in Fig.l. In Fig. 3 the magnitude 
of each arrow gives the absolute value of the transvere component of the total pseudospin 
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sin a and it's angle with the x-axis gives (0^ — (pi). One should note in this regard the 
{(pu — (pi) = A in the other parametrisation Here also one can see the length of the arrow 
does not quite vanish in one of the bimeron centres and makes this construction singular 
at this point. However when we look at the Fig. 4. where the magnitude of each arrow is 
smadownarrow and the angle is again (0„ — (pi) , the magnitude of the arrow vanishes at each 
bimeron centre and the profile is no more singular. This makes Fig. 4 identical to the vector 
plot of the bimeron pseudopsin in the pure layer case [|l5l, |ll3| . 

In Fig 5. we have given a similar vector plot for the spin skyrmion in the upper layer. 
Here the length of the each arrow corresponds to the planar projection of the spin in the 
upper layer (sin6'„) and the its direction gives the azimuthal angle (0„) of the projected 
vector .This picture very clearly points out how the skyrmion winds in the azimutahl plane 
about it's centre at x = b Here also the layer separation d and the starting values of b and 
A are the same as those in Fig. 1 . This set of parameters represents a typical example. 

Lastly, we have evaluated the energy of these solutions for a set of values of the meron 
separation parameter b .The optimal value of b should be obtained by minimising the full 
energy as a function of b. But if we include only the capacitance term and the (pseudo)spin- 
stiffness gradient terms in our energy, these will not lead to a non zero b, i.e the textures 
will want to shrink to zero size. The reason is that under rescaling, the capacitance term 
grows proportional to the square of the scale while the gradient terms are scale invariant. Of 
course a change in b will not result in just an overall rescaling of the solution. The shape of 
the solution will also change since b occurs as a constant multiplying the capacitance term in 
the differential equation (4.5). As a result, the b dependence of the gradient and capacitance 
terms will not be simple, although qualitatively it should still drive the meron separation 
to zero size. This can be seen in Table 1, where we show these energies for ten different 
values of b for a fixed value of other parameters. As expected the sum of both these energy 
contributions decreases strongly with decreasing b, thus driving the texture to zero-size. In 
reality however, these two terms are just the first two terms in the gradient expansion of 
the full energy functional. Higher order terms in the gradient expansion, if included, will 

24 



make our nonlinear differential equation even more difficult to solve, but they can offset this 
tendency to shrink. 

In particular, one prominent higher gradient contribution to the energy is the Coulomb 
interaction between different portions of the topological charge densities. It is given by (see 



where 5p{r) is the CP3 topological charge density given by the integrand of the r.h.s of the 
equation (|3.14| ). 

Inclusion of the contribution of this term into our differential equation for the texture 
will introduce a nonlocal nonlinear term, which will make it very difficult for us to solve 
it numerically. We can however make the following estimate. We can insert our texture 
solution , obtained without the Coulomb term , into the Coulomb energy integral and 
evaluate it as a function of h. This contribution is also shown in the Table 1. As expected, 
the coulomb repulsion energy Ecoui decreases with increasing h. This term would like to 
keep the merons farther apart, and offset the tendency to shrink because of the other terms. 
Thus one may hope to get an optimal bimeron separation at which the sum of all these three 
energy contributions will get minimised as a function of h. 

In Table 1 we have presented our calculation for the layer separation d = .81 where 
different contributions to the total energy are shown along with the their sum, the "total 
energy". For this particular layer separation a distinct minimum is obtained around b = 2.01. 

To see whether this behaviour is common to other layer separations we have plotted in 
fig. 6 the total energy Etotai ^is a fuction of bimeron separation b for a set of layer separations. 
The three sets of points in this figure corespond to three different layer separations, namely 
d equal to .51, .71 and .8/. All these three curves show distict minima for the total energy 
as a function of bimeron separation. We have also provided in Table 2 the size (6) and total 
energy of the optimal bimeron for five values of layer separation. 

Notice from Table 2 that the optimal meron separation b decreases with the increase of 



i 




(4.8) 
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layer separation d . For d = .81 the optimal separation is around b = 2.01 and gradually 
increases to 6 = 2.5/ for d = 0.41 . For the case of pure layer bimeron Brey et. al. also 
found ( a similar behaviour. This behavior could be attributed to the following. By 
lowering the layer separation (decreasing the j ratio) one increases the relative importance of 
Coulomb repulsion among topological charge densities (coming from the intra-layer Coulomb 
energy) to the Capacitance term (coming from the interlayer Coulomb repusion). Hence the 
balancing of the capacitance term by the Coulomb energy will take place at a larger bimeron 
separation, thereby increasing the size of the bimeron. 

We found that resulting energy of these solitons at their optimal sizes varies very little as 
one changes the layer separation. Although bimerons of larger size at lower layer separations 
cost higher capacitance energy , the decrease in the Coulomb energy seems to fully offset 
that. As a result the total energy remains almost the same for different layer separations in 
the range d = 0.4/ to d = 0.8/ that we have studied. 

An important question is whether our spin-pseudospin intertwined solution has a lower 
energy than other candidates among the low-lying excitations. Prominent among these other 
low lying excitations with whom such comparisons have to be done are (i) the particle-hole 
excitations and (ii) purely spin or pseudospin textured solitons. To start with note that in 
the minimal prototype CP3 system (valid in the d = limit ; see |3.8|) the energy is just equal 
to Ecp = 47rp'*(5 (see ref. ( lll|). Now , a pure-spin skyrmion in, say, one of the layers 
can also be written in our CP3 4-spinor notation and will have a CP3 topological number 
Q = 1. So will a bimeron in pseudospin of some spin component. Therefore in the prototype 
CP3 system our spin-pseudospin intertwined soliton with Q = 1 will have the same energy 
as a purely spin or pseudospin textured soliton with Q = 1. The intertwining will not cost 
more energy. All these energies which are equal, are a quarter of that of a particle-hole pair 

— 2 

I ^ 1.25 in units of ^ . However the difference in the 
energies of these various types of topological excitations come from the additional terms in 
the full energy Ec due to capacitance , anisotropy and Coulomb repulsion. As we can see 
from the Table 2 that our intertwined soliton over a range of layer separation has energy 
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around .60 ^. It is encouraging that a pair of these excitations would have somewhat lower 
energy than the particle-hole pair energy of 1.25. 

Of course our computational accuracy is not very high, given that we are limited in 
how many lattice points we can use. One must also improve on the results by solving 
for the texture functions and their energy after including single particle terms due to the 
Zeeman coupling and tunelling . Examples of such calculations can be found in the case of 



u = 2 by Pardes et. al. [|T4| , but not for u = 1 yet to date. Meanwhile, our result for the 
intertwined soliton at z/ = 1 and its energy at best raise hopes that they may be competetive 
as candidates for low lying excitations in double layer systems with spin. 
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TABLE 



Table 1: Different contributions to tfie total energy {Etotai) of spin-pseudospin intertwined 
solitons for a set of 6 at a layer separation of d = 0.8/. Here Egrad refers to the gradient 
energy (isotorpic plus anisotropic), while Ecapa is the capacitance energy and Ecoui is the 
Coulomb interaction energy between topological charge densities . Etotai is the sum of these 

2 

three contributions to the energy. The unit of energy is ^ and the unit of length is I. 



h 


Ecapa 


Egrad 


Ecoul 


Etotai 


4.5 


0.285 


0.261 


0.141 


0.687 


4.0 


0.250 


0.251 


0.152 


0.653 


3.5 


0.229 


0.228 


0.168 


0.625 


3.0 


0.205 


0.223 


0.183 


0.611 


2.5 


0.153 


0.227 


0.217 


0.597 


2.3 


0.143 


0.227 


0.225 


0.595 


2.0 


0.126 


0.203 


0.262 


0.591 


1.8 


0.121 


0.192 


0.290 


0.603 


1.5 


0.104 


0.196 


0.328 


0.628 


1.2 


0.091 


0.200 


0.390 


0.681 


1.0 


0.081 


0.192 


0.456 


0.729 
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TABLE 



Table 2: The size (i.e.the optimal meron separation b) and the total energy (Efotai) at 

2 

different layer separations d. The unit of energy is ^ and the unit of length is I. 



d 


h 


Etotal 


0.8 


2.0 


0.59 


0.7 


2.2 


0.59 


0.6 


2.3 


0.60 


0.5 


2.4 


0.60 


0.4 


2.5 


0.59 
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FIGURES 

FIG. 1. The solution cosq!(77) of the field equations for a set of values for ^.The curves corre- 
spond, as one goes inwards, to (p = O.OOvr, O.SGtt, O.GSvr, O.OOvr respectively with the outermost one 
corresponds to (j) equal to O.OQtt. The layer separation d is equal to 0.6/ and bimeron separation b 
is equal to 2.5/. The value of A in the analytic ansatz is 11 

FIG. 2. The plot of cos a^^o^^^^pj^^) (77) for a set of values for ^.The curves correspond, as one 
goes inwards, to 4> = O.OOvr, 0.367r, 0.637r, 0.907r respectively with the outermost one corresponds to 
(f) equal to O.OQtt. The layer separation d is again equal to 0.6/ and bimeron separation b is equal 
to 2.5/. The value of A in the analytic ansatz is also 1/ 

FIG. 3. This figure gives the magnitude and direction of x-y projection of the total pseudospin 
at different points on the plane. The magnitude of each arrow at a given point is sin a and it's 
angle with the x-axis is (0; — cpu) at that point. The layer separation and the bimeron separation 
are same as in Fig. 1 and 2 

FIG. 4. This figure gives the magnitude and direction of x-y projection of the " pseudospin 
in the down-spin component" at different points on the plane. The magnitude of each arrow at a 
given point is sinaj^ and it's angle with the x-axis is at that point. The layer separation and 
the bimeron separation are same as in Fig. 1 and 2 

FIG. 5. This figure gives the magnitude and direction of x-y projection of the spin in the upper 
layer at different points on the plane. At each point the magnitude of the arrow gives sin 0^ and 
it's direction with x-axis gives 0„ at that point. The layer separation and the bimeron separation 
and initial A are the same as in the earlier figures 

FIG. 6. This figure gives a plot of the total energy E(total) as a function of the bimeron 
separation b for three different layer separations, namely d = 0.5/, 0.7/ and 0.8/ The unit of energy 
IS agam 
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